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ABSTRACT 

This paper describes an imaging microscopic technique based on heterodyne digital holography where subwavelength- 
sized gold colloids can be imaged in cell environment. Surface cellular receptors of 3T3 mouse fibroblasts are 
labeled with 40 nm gold nanoparticles, and the biological specimen is imaged in a total internal reflection con- 
figuration with holographic microscopy. Due to a higher scattering efficiency of the gold nanoparticles versus 
that of cellular structures, accurate localization of a gold marker is obtained within a 3D mapping of the entire 
sample's scattered field, with a lateral precision of 5 nm and 100 nm in the x,y and in the z directions respectively, 
demonstrating the ability of holographic microscopy to locate nanoparticles in living cells environments. 
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1. INTRODUCTION 

With recent developments in the fields of nanotechnology and modern optical microscopy, the use of nanometric 
particles as biomarkers in biological specimens has been rapidly increased. Colloidal gold nanoparticles have 
gained popularity over other biomarkers for several reasons. First, gold nanoparticles provide high scattering 
efficiencies. 1 Second, unlike fluorescent markers, they are immune to photo bleaching, and they are potentially 
non-cytotoxic, which provides a great advantage over semiconductor nanoparticles. 2 Because of these properties, 
the potential for the use of gold nano particles as biomarkers for live cell imaging using photothermal tracking 3-5 
and light scattering spectroscopy 6 has been demonstrated. Detection of live oral cancer cells using gold nanopar- 
ticles have also been obtained. 7 In this paper, we show the potential of using Digital Holographic Microscopy 
(DHM) as a powerful tool to detect and localize, in three dimensions, gold nanoparticles attached to living cells. 

In digital holography a charged coupled device detector (CCD) records the interference diagram, called 
hologram, of two waves: the reference wave and the object wave. 8 The interference pattern is then numerically 
reconstructed to obtain information about the object wave field. A single hologram can be used to reconstruct 
the wave field at virtually any plane, thus allowing one to obtain both amplitude and phase information of 
the whole object. Over the past few decades, different configurations and techniques have been introduced in 
digital holography, resulting in precise phase shift, 9 high resolution, 10-12 and high sensitivity measurements. 13 
Some of these techniques involve improvement of the experimental setup 14 and data acquisition methods, while 
others involve additional numerical treatments of the acquired holograms and modification of the numerical 
reconstruction procedure. 15 In this vision, Leith and Upatnicks 16 suggested an off-axis configuration holography, 
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Figure 1. Schematic representation of the coupling between 3T3 cells and gold beads 

in which the combination of the two beams is done after angularly tilting the propagation direction of one of 
the beams with respect to the other. Phase-shifting interferometry adapted to digital holography 17 has been 
widely used, and our group has proven the potential efficiency of combining off-axis and phase shifting with a 
spatial filtering technique for reducing noise, 13 and eliminating images aliases and overlapping. 18 The progress 
in various digital holographic techniques has dramatically increased its applications to a large variety of fields, 
and especially in the fields of cell biology and biomedical microscopy. 19-22 

Holography has now proved its ability to localize scattering nanoparticlcs in 3D, as shown by Atlan et al. 
in a recent paper, 23 either for fixed particles spin coated on a glass substrate or in free motion within a water 
suspension. More recently, Absil et al. 24 have shown that heterodyne holography also allows the photothermal 
detection of small nanoparticles. 

However, in order to apply these techniques to biological specimens, important issues must be considered. 
In biological samples, the particle holographic signal is superimposed with the light scattered by cell refractive 
index fluctuations, which yield a speckle field. While in our case this speckle is a parasitic signal, in many 
other situations, like in Dark Field microscopy, or in Differential Interference Contrast (DIC) microscopy, 25 this 
speckle is the main source of contrast used to image the cell. In this paper, we study the possibility of 3D 
holographic imaging in a biological context and we image for the first time 40 nm gold nanoparticles attached 
to living cells using DHM. Since the cell scattered speckle field cannot be avoided, it is important for future cell 
labeling applications, to scale the particle signal with respect to the scattered speckle. We show here that the 
amplitude of the 40 nm gold particle signal we get is much larger than the cell scattered field. 

2. SAMPLE PREPARATION 

We use NIH 3T3 mouse fibroblasts (quoted as 3T3 cells in this paper) with integrin surface receptors labeled with 
40 nm gold particles. Streptavidin-coated gold nanoparticles were attached to the cell surface integrin receptors 
via biotin and fibronectin proteins: see Fig. [T] Streptavidin and biotin are very well known for their strong 
affinity towards each other, and fibronectin, an extracellular matrix protein, has the property of interacting 
specifically with cell surface receptors of integrin family. 

Fibronectin proteins (fibronectin from bovine plasma, Sigma, St Louis, MO) were labeled with biotin by using 
EZ-Link®Sulfo -NHS-LC-Biotin according to the provider protocol (Pierce, Rockford, IL). The final concentra- 
tion of biotinylated-fibronectin solution was 0.447 mg/mL. We used a solution of gold nanoparticles of 40 nm 
average diameter pre-coated with streptavidin. The streptavidin/gold conjugates (BioAssay, Gentaur, France) 
were rinsed twice with 1 M PBS (Phosphate Buffered Saline) (pH = 7.25) and 10 fiL of the gold solution was 
added to 990 fiL of the same PBS buffer to form 1 L dilute solution. Then the dilute solution was incubated with 
50 /iL biotinylated-fibronectin solution for four hours at room temperature to allow streptavidin-biotin bonding. 
The functionalized beads solution was stored at 4°C and used within 24 hours after preparation in order to 
ensure maximum functionality. Before every use, the functionalized beads solution was sonicated. 
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Figure 2. Experimental setup. AOM1, AOM2: acousto-optic modulators; M: mirror; MO: microscope objective NA = 
0.5; A/2: half wave plate; BS: beam splitter; PBS: polarizing beam splitter; CCD: CCD camera; Er\ reference wave; Eo- 
object illumination wave; E: scattered wave; 8: angular tilt; NDi, ND2: neutral density filters; z = 0: CCD plane; z = zo: 
CCD conjugate plane with respect to MO. 



3T3 cells were cultured in Duclbecco's modified Eagle's medium (DMEM Gibco, Invitrogen, Carlsbad, CA) 
supplemented with 10% fetal calf serum (PAA Laboratories GmbH) on 32 mm diameter glass cover slips at 37°C in 
a 5% C02 atmosphere, 48 hours before the observation. The cover slips were coated with fibronectin (fibronectin 
from bovine plasma, Sigma, St Louis, MO) in advance for optimum cell growth. After being incubated for 24 
hours, a monolayer of 3T3 fibroblasts was immersed into a solution composed of DMEM (2 mL) plus 500 /iL 
of the functionalized beads solution. The integrin-fibronectin link is created at this level allowing the cells to 
attach, on their surface, the functionalised gold nanoparticles. 

The coverslip containing adherent 3T3 cells tagged with 40 nm gold nanoparticles was mounted on a specific 
observation chamber. During the observation, cells were kept in DMEM F12 medium with 0.5% fetal calf serum 
with no phenolred, to fulfill optimum survival conditions. The efficiency of biotinylation protocol is verified by 
measuring the level of biotin incorporation on an HABA [2-(4' -Hydroxyazobenzene) Benzoic Acid)] quantitation 
assay. Average number obtained of biotin molecules per fibronectin is 2.5. 

3. EXPERIMENTAL SETUP 

A schematic representation of the experimental setup is shown in Fig. [2j The illumination source is a single mode 
near infrared laser diode emitting at A=785nm (50 mW, 90 mA current). A polarizing beam splitter cube (PBS) 
is used to split the original illumination laser light into two beams, a reference beam (complex field Eji, frequency 
f R ) and an object illumination beam (complex field Eo, frequency f ) forming the two arms of a Mach-Zchndcr 
interferometer. A combination of a half wave plate (A/2) and two neutral density filters (NDi, ND2) is used to 
prevent the saturation of the detector by controlling the optical power traveling in each arm. Two acousto-optic 
modulators (AOM1, AOM2) driven around 80 MHz with a selection of the first order of diffraction, shift both 
frequencies at respectively Jaomi and Jaom2 ■ 

The object beam illuminates the sample, prepared as described in section [2j by provoking total internal 
reflection (TIR) at the cell/bead-air interface. The evanescent wave locally frustrated by the beads and cells 
gives off a propagating scattered wave (complex field E), which is collected by a microscopic objective (MO, 
50 x magnification, NA=0.5, air). A beam splitter is then used to combine the scattered object wave and the 
reference wave which is slightly angularly tilted (6 ~ 1°) with respect to the propagation axis of the object wave 
in an off-axis configuration. The half wave plate (A/2) on the object illumination arm aligns the polarization of 



the corresponding beam allowing an optimal interference of the beams. A CCD camera (Roper Cascade 512F, 
512x512 square pixels of 16/im size, exposure time 100 ms, frame rate fccD — 8 Hz) detects the interference 
pattern (hologram) and sends it to a computer. The hologram is then numerically treated and the complex field 
E(x,y,z) is reconstructed numerically. 

4. HOLOGRAPHIC ACQUISITION AND RECONSTRUCTION 

A four-phase demodulation method is used to record holograms. This method consists in acquiring a sequence 
of images with a relative phase shift Aip = tt/2 between two consecutive frames. To get an accurate phase shift, 9 
the reference wave is frequency shifted by tuning the two acousto-optic modulators AOM1 and AOM2, 14 so that 
the heterodyne beat frequency is: 

fcCD 

A/ = fAOMl — JAOM2 — — t — (1) 

where fee d =8Hz is the frame rate frequency of the CCD camera. The camera records a sequence of 32 frames 
Io,...,l3i, and the object field E on the CCD plane (z = 0) is given by: 

M 

E(x,y,z = 0) = J2j ni n (2) 

71=1 

where j 2 = —1, and M is the number of frames used for the reconstruction, which is equal either to 1 or to 32 in 
the experiment we present here. In Eq{2] the coordinates x, y (with < x, y < 511) are integers, which represent 
the pixel location within the CCD plane. The pixel size is then the physical CCD pixel size, i.e., 16 /im. 

Since we image the sample through a microscope objective (MO), the measured hologram represents the field 
E{x,y,z — zq) within the CCD conjugate plane z = zq with respect to MO, i.e., the plane, whose image is on 
focus on the CCD detector. In that case, one must compensate the phase curvature, the phase tilt and the 
enlargement factor that are related to the presence of MO. 26 We have thus: 

M 

E(x, y, z = z ) = e m^+K v y) e M(* W) J- j n I n (3) 

n=0 

where (K Xl K y ) and A are the tilt and lens parameters respectively. These parameters are determined by 
reconstructing the image of the microscope objective output pupil using the one Fourier transform method. 27 
The lens parameter is close to the lens parameter that is used in the pupil reconstruction, and the tilt parameter 
to the translation that pushes the pupil in the centre of the reconstructed image. The magnification factor 
of the conjugate plane is measured by imaging a USAF target. We get, in the conjugate plane, pixel sizes of 
Sx, Sy = 177 nm. 

The object field E(x,y,z) is calculated in the vicinity of conjugate plane (i.e. for z ~ zq) by the angular 
spectrum method, which requires two Fourier transforms. 10 ' 14,28 This method is chosen in order to keep the 
pixel size constant in the reconstruction. The reconstruction is done for 512 different reconstruction distances 
z = zo + mSz where Sz = 177 nm and m = —256... + 255. By this way, we get 3D volume images with 
512 x 512 x 512 voxels, with the same pixel size (177 nm) in the 3 directions x, y and z. 

5. EXPERIMENTAL RESULTS 

We have considered a 3T3 cell with intcgrin cellular surface receptor tagged with a single 40 nm gold particle. 
Since we are interested in tracking the particle, the camera exposure time is an important issue, and we will first 
reconstruct the 3D image of the sample using one CCD frame. Then, to improve the Signal to Noise Ratio and 
to confirm our results, we will reconstruct the 3D image of the sample from the whole sequence of 32 frames. 
Figure [3] (a) shows a direct image of the sample under white light illumination. While the cell is well defined on 
the left hand side of the image, the attached gold particle cannot be seen. 

Figure [3] (b) shows the reconstructed intensity image of the sample in a colored logarithmic scale. The 
reconstruction is done using 1 CCD frame (M = 1) with an exposure time of 100 ms. The reconstruction plane 




Figure 3. Fibroblast cell tagged with a 40 nm gold particle, (a) Direct image under white light illumination, (b) 
Reconstructed Holographic intensity I image. The 40 nm gold particle is marked with a white arrow. The color scale 
corresponds to 6 < ln(J) < 15). (c) Volume view of the 3D reconstructed data (512 x 512 x 512 voxels; voxel size 177 nm 
in all directions). The holographic reconstruction is made in (b) and (c) from 1 CCD frame with an exposure time of 100 
ms. 
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Figure 4. A fibroblast cell tagged with a 40 nm gold particle using 1 frame acquisition, (a) 3D linear surface plot of the 
intensity image reconstructed from a single frame, (b) Experimental linear plot of cuts made within the intensity signal / 
along the three axis x,y and z at the brightest voxel of Fig[3](b) corresponding to the gold nanoparticle location. Curves 
1, 2 and 3 corresponds to the x, y and z axis respectively. The horizontal black dashed line represents the half maximum 
of the curves. Curves marked with arrows are zooms of the corresponding curves from maxima to half maxima. Light 
grey curves correspond to the different ideal theoretical curves. 



is z = zq + 12.2 /im i.e. n z = 325 where n z is the z index on the 3D matrix, which is the plane where the brightest 
point is detected. It also corresponds to the cell plane. The background signal is the red and blue triangular 
structure, with a few speckle yellow high spots. It is due to the light scattered by the inhomogeneities of the 
cell refractive index, and thus corresponds to the cell imaged in Fig{3] (a). The gold particle corresponds to the 
brightest spot in Fig|3] (b). It is located on the left bottom of the cell, and marked with a white arrow. 

To illustrate the ability of our technique to localize the gold particle in 3D, we have displayed in Fig. [3] (c) 
a volume view of the 512 x 512 x 512 reconstructed intensity image 3D data, i.e., we have plotted I(x,y,z) — 
\E(x, y, z)\ 2 for different reconstruction z distances. The bright light bluish zone on the 3D image corresponds 
to the 40 nm gold particle signal. This wavefield has the shape of a cigar roughly oriented along the microscope 
objective axis (z axis). 

To better visualize the 40 nm gold particle, we have displayed, on Fig. [4] (a), a 3D linear surface plot of 
the corresponding logarithmic scale colored image ofFig|3](b). As seen, the gold particle signal is much higher 
than the background signal that corresponds to the speckle scattered by the cell. The particle can thus be easily 
located. 

To perform a quantitative analysis of the precision we expect in the 3D localization of the gold particle, we 
have considered the 3D volume intensity image of Fig(3](c), and we have made cuts along the 3 axis (x, y and z) 
at the gold particle crossing point (i.e. where the three axis cross the highest intensity voxel of the 3D image). 
The intensity signals I along the 3 axis are plotted on Fig|I](b), curves 1 (red), 2 (blue) and 3 (green) correspond 
to x, y and z axis respectively. The horizontal black dashed line indicates the half maximum of the curves. In 
order to measure the curves widths, we have displayed zooms of the curves from maximum to half maximum. 
The individual pixels are visible on curves 1 and 2 zooms (x and y axis). The Full Width at Half Maximum 
(FWHM) is about 6 pixels in the x and y directions (6 x 177 = 1060 nm), and about 60 pixels in the z direction 
(10.6 /xm). 

To compare the resolution obtained here with that expected from the NA of the microscope objective, we have 
computed the wavefield Ej(x, y, z) that is expected in the ideal case for a NA=0.5 objective, and we have made 
cuts along x, y and z directions of the 512 x 512 x 512 ideal cube of data for \Ei(x, y, z) | 2 . The corresponding 
theoretical cuts are plotted with the zooms in heavy grey line. Like for the experimental cuts, the horizontal axis 
is in pixels. In the x and y directions, the obtained FWHM is close to the expected one. Yet, in the z direction 
the shape of the curve is not symmetric, and the measured FWHM is about 3 times larger than the expected 
FWHM in the ideal case. This discrepancy may be related to the parasitic speckle background signal. 





Figure 5. A fibroblast cell tagged with a 40 nm gold particle using 32 frames acquisition, (a) 3D linear surface plot of 
the 32 frames reconstructed holographic intensity image, (b) Reconstructed Holographic intensity 7 image. The 40 nm 
gold particle is marked with a white arrow. The color scale corresponds to 6 < ln(7) < 15). 



Using a parabolic approximation for the local field, the location of the gold particle can be calculated by 
fitting the data points that are above half maximum. The accuracy of the measurement made by this method 
is ±5 nm in the x and y directions (uncertainty given by the fit software Gnuplot). It corresponds also to the 
variation of the particle when the number of data points used for the fit is increased from 8 to 7. In the z 
direction, the measured curve is not perfectly symmetric, and so the fit technique is not very accurate because 
the result of the fit strongly depends on the data points that are used in the calculation. 

In order to visualize the gold particle with a better Signal to Noise Ratio (SNR), we have also reconstructed 
the holographic images using a sequence of M — 32 acquired CCD frames. In this case, the total exposure time 



(a) shows the 3D linear surface 
5] (a) is similar to Fig|3] (a) but 



is 32 x 100 ms=3.2 s, while the measurement time is 32/fccD = 4 s. Figure p 
plot of the reconstructed intensity image in the z = z + 6.1/im plane. Figure 

obtained with 32 frames, instead of 1 frame. Here again, not only the gold particle signal is much higher than 
the background signal, but the ratio of the particle signal versus background signal (light scattered by the cell) is 
increased with respect to the ratio obtained for Figj4](a), and the visibility of the gold particle is improved. This 
visibility improvement is confirmed by Fig. |5](b), which shows the reconstructed intensity image of the sample 
in a colored logarithmic scale. The cell signal, corresponding to the triangular bluish and reddish structure, is 
still visible. The particle signal in yellow is visible too, but with a much better contrast than in Figj4](a). 
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Figure 6. (a) Experimental linear plot of cuts made within the intensity signal I along the three axis x, y and z at the 
brightest voxel of Fig[5] (b) corresponding to the gold nanoparticle location. Curves 1, 2 and 3 corresponds to the x,y 
and z axis respectively. The horizontal black dashed line represents the half maximum of the curves. Curves marked 
with arrows are zooms of the corresponding curves from maximum to half maximum. Light grey curves correspond to the 
different ideal theoretical curves, (b) Intensity plot of the cut along the x axis in logarithmic scale 1 at the nanoparticle 
location and 4 without illuminating the sample 



Figure [6] (a) shows the cuts made along the 3 axes (x, y and z) obtained for the 32 frames data. Here 
again, we have displayed zooms of the three cross sections from maxima to half maxima as well as zooms of 
the theoretical curves. These cross sections are similar to the cross sections obtained in the case of 1 frame 
acquisition, but in this case (32 frames) the heavy grey lines corresponding to the theoretical zoomed curves are 
closer to the experimental curves. This confirms our hypothesis of the effect of the background parasitic signal 
(which is much lower here, as seen by comparing Fig{5](a) with Figj4](a)). Here again, the location of the particle 
can be calculated by fitting the measured data point with a parabola. The accuracy is still ±5 nm is the x and 
y directions, but since the curve 3 is now roughly symmetric, the fit can also be performed in the z direction, 
and the particule can be located with an accuracy estimated to ±100 nm in the z direction. 

To better compare the particle signal with respect to the light scattered by the cell (background speckle), 
we plotted the intensity signal along x (which corresponds to curve 1 of Fig{6] (a) ) in logarithmic scale on Figj6] 
(b). As seen, the particle signal is 2 to 3 orders of magnitude larger than the background signal. It has the same 
order of magnitude that the signal obtained in a control experiment of 40 nm gold particles within a water and 
agarose suspension. We plotted also (black curve 4) the background signal obtained without illuminating the 
sample. This background corresponds to the ultimate sensitivity limit of the holographic detection, 29 i.e. to the 
shot noise on the local oscillator beam. As seen, the shot noise background is about 100 x lower than the cell 
background. This means that if the exposure time is reduced by the same factor (from 3.2 s down to 32 ms), one 
expects to get roughly the same Signal to Noise Ratio (SNR~ 100) on the particle signal. Since the illumination 
area (about 1 mm 2 ) is much larger than the imaged region, it is possible, by focusing the illumination, to reduce 
further the exposure time without significally decreasing the SNR. 

6. CONCLUSION 

In conclusion, we have reported in this article a digital holographic microscopy technique capable of imaging 40 
nm gold nanoparticles attached to living cells surface receptors. Results show that this method can be effectively 
used to distinguish nanoparticles from the cellular structures. We show in particular that the 40 nm gold particle 
scattered signal is much larger in intensity than the field scattered by the cell structures. While a single acquired 
image is sufficient to distinguish the signal of the particle from the cell scattered field, better signal to noise ratio 
and z-positions estimations can be obtained if the reconstruction is done by averaging on a sequence of several 
acquired images. Overall, we showed that we can use our technique for the 3D mapping of the sample's full 
scattered field using a single recorded hologram and without any mechanical scanning. 
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